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ABSTRACT 

During the initial data reduction of the Wisconsin H- Alpha Mapper (WHAM) Ha Sky 
Survey, we have discovered several very long (~ 30°-80°) filaments superimposed on the diffuse 
Ha background. These features have no clear correspondence to the other phases of the 
interstellar medium revealed by 21 cm, X-ray, IR, or radio continuum surveys, and they have 
no readily identifiable origin or source of ionization. In this letter, the data for two of these 
faint (Iua ~ 0.5-1.5 R) structures are presented. The first is an 80°Tong, 2°-wide arch that 
extends nearly perpendicular to the Galactic plane at I = 225° and attains a maximum latitude 
of +51° near I = 240° before reaching the southern boundary of our survey map at I = 270°, 
b = +42°. The vertical portion of this feature between b = +10° and +25° is associated with 
a single radial velocity component centered at wlsr = +16 km s" 1 with a full width at half 
maximum of 27 km s" 1 . A decrease in the velocity is observed from b = +33° through +48° 
as the feature arches toward higher Galactic longitudes. At this end, the emission component 
is centered near wlsr = ~20 km s" 1 . Where this feature appears to meet the Galactic plane 
near I = 225°, it is directly above the H II region surrounding CMa R1/OB1. A second filament 
consists of a - 25°-30°-long arc spanning £ = 210°-240° at b = +30° to 40°. The radial velocity 
of this feature increases systematically from km s" 1 at I = 215°, b = +38° to +18 km s _1 at 
£ — 236°, b — +28°. Both features have rather constant intensities along their entire lengths, 
ranging from 0.5-1.5 R (EM = 1-3 cm" 6 pc) with no obvious trends. 

Subject headings: ISM: structure — Galaxy: structure — Galaxy: halo 



1. Introduction 

The Wisconsin H-Alpha MapperQ (WHAM) survey is providing the first velocity-resolved map of the 
Ha emission from our Galaxy's diffuse interstellar medium. The combination of WHAM's sensitivity (< 0.1 
R) and velocity resolution (12 km s" 1 ) reveals new details about the large-scale structure and kinematics of 
the Warm Ionized Medium (WIM). Studies of faint emission structures in the WIM may help us understand 
the distribution and morphology of the ionized gas and may even provide direct evidence of ionization 
sources. In this letter, we present the discovery of two large, faint filaments in our Galaxy. 
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2. Observations 

WHAM consists of a dedicated 0.6-m, all-sky siderostat connected to a 15-cm, double-etalon 
Fabry-Perot spectrometer. The optical design delivers a spectrum covering a 200 km s _1 radial velocity 
interval from a one-degree circular patch of sky. For the Ha sky survey, approximately 35,000 spectra were 
obtained above 5 = —20° with the one-degree diameter beam centered on a 0?98 x 0?85 grid in £ and b. The 
integration time for each one-degree pixel was 30 seconds, resulting in a signal-to-noise ratio of about 30 in 
background continuum regions of the spectrum. This allows the detection of sources as faint as 0.1 R (1 R 
= 10 6 /4-7r ph cm -2 s _1 ster -1 or 2.4 x 10~ 7 ergs cm -2 s _1 ster -1 at Ha) in a single exposure. The etalons 
are configured to provide a spectral resolution of 12 km s -1 , which is sufficient to identify and remove 
atmospheric lines in the spectra and fully resolve the Ha emission-line profiles from warm gas (T ~ 8000 
K) in the Galaxy. A detailed description of the instrument design is given by Tufte (1997). 

Presented here are nearly 3800 of these spectra extracted from the WHAM Ha Sky Survey. They cover 
the region of the sky from I = 200° to 270° and b = -10° to +60° above 5 = -20°. Most of these spectra 
were observed by WHAM from 1997 February through 1997 May. These observations are designated "01" 
below. A small portion (145 spectra) from I = 253° to 270° and b = +47° to +60° was observed in January 
1998 and are labeled "02" below. Standard WHAM data reduction, including bias subtraction, reflected 



ring subtraction, annular-summing, and flat fielding, was used to create these spectra (|Haffncr et al. 1998| ) 



Atmospheric emission lines are a significant source of contamination in WHAM spectra. The brightest 
of these, with intensities ranging from 2-13 R, is the geocoronal Ha line emitted in the earth's upper 
atmosphere. Since the earth's orbital velocity changes the location of the velocity frame of the Local 
Standard of Rest (LSR) with respect to the geocentric velocity frame, most directions in the sky were 
observed when the geocoronal line was separated by at least 20 km s _1 from the LSR. This separation 
makes it easy to resolve the relatively narrow (6-8 km s -1 FWHM) geocoronal line from the Galactic Ha 
emission and minimize its effect. The data presented here have velocity separations of the geocoronal line 
from the LSR that range from -24 km s" 1 to -42 km s" 1 (01) and +14 km s" 1 to +25 km s _1 (02). A 
single Gaussian component has been fitted and subtracted from these spectra to remove the geocoronal line. 

Two other much fainter atmospheric emission lines, believed to be extremely weak OH emission, are 
present in the spectra. The first is located 71 km s _1 to the red of the geocoronal line («lsr = +27 to 
+45 km s _1 in 01); it has a width of 23.5 km s _1 and varies in intensity from 0.05 to 0.36 R. Deeper 
observations and the wide line width suggest that this feature may be a pair of atmospheric lines. The 
second is located 38 km s _1 to the blue of the geocoronal line (i>lsr = —64 to —82 km s -1 in 01); it has a 
width of 9.4 km s _1 and varies in intensity from 0.04 to 0.13 R. For this study, the red sky line provides the 
most contamination near the LSR in 01, and it has been removed by approximating it as a single Gaussian 
component. Variations in the line intensities are most severe from night to night. Our corrections account 
for intensity variations over 30-minute periods, which are typically 0.02 R. Other faint atmospheric lines are 
present in 02 even further to the blue from the geocoronal line. These features will not be discussed here 
since they do not affect our study near the LSR. 

A first-order polynomial fits the sky background well within our 200 km s _1 bandpass. The typical 
level of the subtracted background is approximately 0.03-0.07 R (km s -1 ) -1 or 1.4-3.2 R A" 1 , which is 
dominated by atmospheric continuum emission above \b\ = 30°. Individual stars within the one-degree 
beam brighter than sixth magnitude contribute significantly to the background in several directions. Some 
of the spectra in these directions also contain an Ha absorption line from the stars. These pointings are 
most noticeable in Figure [l] as isolated one-degree diameter depressions in the otherwise smooth emission 
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at higher Galactic latitudes. 

Absolute intensity calibration is accomplished by frequent observations of nebular sources. All WHAM 
intensities are tied to the absolute intensity measurement of the North American Nebula (NAN; NGC 
7000), 850 R within a 50' beam centered at £ = 85?60, b = -0°71, as determined by Scherb (1981). Since 
we have yet to apply the gamut of intensity calibrations to our data, relative intensities quoted here are 
accurate to about 10%. The systematic uncertainty in the absolute scale (i.e., the intensity of NGC 7000) 
is believed to be approximately 15%. 



3. Results 

Figure [l] displays four velocity integrated images of Ha emission from our selected region of the sky. 
The two filaments of interest are most easily seen in the narrow (12 km s _1 ) velocity interval images, 
although their overall structure can also be seen on the "broad-band" map (Figure |l|a). The vertical portion 
of the longest filament (hereafter Filament 1) is most prominent in the map centered on ulsr = +20 km s _1 
(Figure |l]d), stretching nearly perpendicular to the Galactic plane along £ — 225°. It then extends toward 
higher Galactic longitudes, turning over near b — +50° while shifting to more negative velocities in Figures 
[j]c and [l]b. A second long filament (Filament 2) can been seen as an arc extending across the images near 
b = +30° to +40°. Its peak emission moves from £ — 220° in the i>lsr = — 20 km s _1 map (Figure |l|b) 
to £ = 235° in the t+SR = +20 km s _1 map (Figure |]d). Other filaments are also apparent, including the 
feature running from £ = 205°, b = +25° to £ = 220°, b = +30°. In this letter we focus primarily on the 
80°-long Filament 1. 

An £-v map is presented in Figure || along b = +25°. From this representation, Filament 1 appears to 
be a narrow (<; 2°) feature at £ = 225° that is offset in velocity from the smooth emission present near the 
LSR at other Galactic longitudes. An analysis of the spectra reveal that this velocity offset is about +18 
km s _1 . The spectrum of Filament 1 averaged over six of the WHAM pointings taken around £ = 225?5, 
b = +25?5 is displayed in Figure || as a dashed line. Because the velocity separation of the two components 
(Aw ps 20 km s" 1 ) is less than the typical Doppler width of WIM Ha lines (FWHM ps 20 - 30 km s" 1 ), a 
straightforward two-component fit of this spectrum is not well constrained. Instead, we attempt to isolate 
the filament from gas near the LSR by subtracting an estimate of the general WIM contribution in this 
direction. The spectra of two "background" regions taken at the same Galactic latitude are also displayed 
in Figure || as dotted lines. These directions are averages of nine pointings around £ = 232?0, b = +25?5 
and 12 pointings around I = 220?5, b — +25?5. An average of these spectra is subtracted from the average 
filament spectrum to produce the final filament-only spectrum in Figure [| (solid line). The resulting 
component is fit well by a single Gaussian with a mean of +18 ± 2 km s _1 , a FWHM of 27 ± 6 km s _1 , and 
an intensity of 1.2 ± 0.2 R. 

Using this method, we have determined fit parameters for several directions along Filament 1 and 2 and 
present the results in Table [1} Regions where the filaments cross are deliberately excluded. The accuracy 
of the filament component parameters depends primarily on how well the background regions represent the 
true emission near ulsr — km s" 1 in the filament direction. The errors of the parameters listed in Table |l| 
are conservatively estimated by examining subtractions of different background directions from the filament 
spectra. The fitting errors are an order of magnitude smaller and do not contribute significantly. Note the 
velocity gradients in both filaments with increasing Galactic longitude, decreasing from ^lsr +18 km s _1 
to —25 km s _1 in Filament 1 and increasing from ulsr km s _1 to +18 km s _1 in Filament 2. Also note 
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the nearly constant Ha intensities and line widths. 

4. Discussion 

The intensities in Table [l] are converted to emission measures using the formula: 

EM (cm" 6 pc) = 2.75 T 4 a9 J Hq (R) e 2 - 2E ^- y \ (1) 
where T4 represents the temperature of the emitting gas in units of 10 4 K. With typical values for the 



temperature of the WIM gas, T4 = 0.8 ( Reynolds 1985 ), an Ha intensity of 1 R corresponds to an EM of 
2.25 cm -6 pc. Since we are examining mostly high-latitude features, we expect corrections for interstellar 
extinction to be small. No such correction has been applied to the data presented here. 

To derive further physical parameters for the filaments, we must make an assumption about the 
distances to them. Some evidence suggests that the lower part of Filament 1 is at a distance of about 1 
kpc, based on its radial velocity and possible association with CMa R1/OB1. The lower portion has a 
radial velocity of +16 km s -1 , which is consistent with a kinematic distance of approximately 1 kpc at 



this Galactic longitude for a value of Oort's Constant A = 16 km s -1 kpc" 1 (Mihalas & Binney 1981). 
Unfortunately, this argument is weakened by the radial velocities of the filament's high-latitude portion, 
which are "forbidden" in the Galactic rotation model. However, in Figure |l|d («lsr = +20 km s -1 ), the 
vertical portion of Filament 1 appears to end in the prominent H II region surrounding the star-forming 
association CMa R1/OB1. In fact, the center of the brighter, arc-shaped, Ha emission features in this 



H II region, as seen on the Palomar Sky Survey red plates, is at £ = 225°, b = —1° ( [Reynolds & Ogden 
1975). Also, the velocity centroid of the H II region (Table g) is similar to that of the first 15°-20° of the 
filament away from the plane. Unless the agreement between the spatial and velocity information of these 
two regions is a coincidence, we can use the distance of the H II region as an estimate of the distance to the 
filament. The distance to the CMa OBI association has been determined photometrically by Claria (1974) 
to be 1150 ± 140 pc, consistent with a kinematic distance of 800-1100 pc for the Ha emitting gas in the 



region ( [Reynolds fc Ogden 1978|) . In this letter, we adopt a distance of 1 kpc for Filament 1 but include 



d\ = d/l kpc in the derived parameters. 

At a distance of 1 kpc, the vertical extent of Filament 1 (t « +51°) translates to a vertical height 
above the Galactic mid-plane, Z, of 1200 d\ pc and a width (« 2°) of 35 d\ pc. If we assume the filament is 
a cylinder of uniform density gas of width L, then the emission measure, 

EM (cm -6 pc) = I n\dl = n 2 e L, (2) 
Jo 

of the feature can be used to estimate the density. For Filament 1, where L = 35 d\ pc and 
EM ss 1.1 cm~ 6 pc, n e — 0.18 d 1 l ' 2 cm~ 3 . The typical column density through the filament is 
N e — 1.9 x 10 19 d 1 ^ 2 cm" 2 . The mass of the material in the filament's vertical section can be estimated 
by 1.4 run Z L 2 n e = 9.2 x 10 3 d^ 2 M Q , where the factor 1.4 is a correction for helium. If the filament 
is photoionized, our observed EM implies an incident ionizing flux of at least as EM = 1.0 x 10 6 ph 
cm -2 s -1 . The hydrogen recombination rate within the filament is given by as n 2 — 1.0 x 10 — 14 d^ 1 
cm~ 3 s -1 , implying that the power required to sustain this rate throughout the volume of the filament is 
(a B n 2 Z L 2 ) x 13.6 eV = 9.4 x 10 36 d 2 erg s" 1 . The value of a B used in these calculations, 3.10 x 10~ 13 
cm 3 sec -1 , is interpolated from Osterbrock (1989) assuming a gas temperature of 8000 K. 
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The narrow shape and location of Filament 1, particularly the coincidence in location and radial 
velocity of one end of the filament with an energetic source in the plane (CMa R1/OB1), suggest one 
possibility for its origin — a jet-like ejection from the association. However, given the length of the filament 
at the distance of CMa R1/OB1, the apparently low velocities associated with it (Table [l]), and the near 
constant Ha intensity profile along its entire length, it is difficult to reconcile a scenario in which ionized 
hydrogen is being ejected from CMa R1/OB1 with the gas's short recombination times. A parcel of ionized 
gas at these temperatures and densities will recombine in t r = (ag «e) _1 = 1.8 x 10 13 d\^ 2 s = 5.7 x 10 5 d\ 
yr. To reach the observed height and remain ionized, the gas would need to be ejected at speeds in excess 
of Z/t r = 1400 d\ /2 km s _1 . We see no evidence for such speeds in our data. Furthermore, we would expect 
a significant gradient in the filament's intensity as a function of height above the Galactic plane, since gas 
at larger distances from the plane has had more time to recombine. 

If another source is responsible for maintaining the ionization of Filament 1, an ejection scenario could 
still be plausible. If the filament's arc shape is caused by free-fall of ejected gas back to the Galactic plane, 
an initial velocity of about 70 km s _1 is required to reach the observed height above the plane at the 
distance of CMa R1/OB1. Such speeds agree better with the measured filament component velocities in 
Table particularly when projection effects are considered. However, the filament's length would then 
require that the ejecting source be more than 3 x 10 7 yr old. This number can be compared to estimates of 
3 x 10 6 yr for the age of CMa OBI by Claria (1974) and 7.6 x 10 5 yr for the age of a supernova put forth by 
Herbst & Assousa (1977) as being responsible for generating the ring-shaped Ha nebulosity in the vicinity 
of CMa OB1/R1 and for being the progenitor of star formation in the Rl association. In this scenario, a 
diffuse ionization source is probably required to maintain a constant intensity along the filament, since the 
growing distance from CMa R1/OB1 would produce an intensity gradient if it were the sole source. The 
efficiency of leaking Lyman continuum radiation from the disk is a longstanding problem for ionizing a 
thicker layer of the Galaxy. However, Bland-Hawthorn & Maloney (1998), Dove & Shull (1994), and Miller 
& Cox (1993) have argued that a substantial ionizing flux can diffuse into the Galaxy's halo from the disk. 
The required ionizing flux of at least 1.0 x 10 6 ph cm -2 s _1 is consistent with these models. 

Several other phenomena may produce large, faint filaments. Large-scale, kinematic Galactic structures 
such as superbubbles, chimneys, or worms have been popular explanations for filamentary structures seen 



in H I as well as Ha {e.g. Koo, Hciles, & Reach 1992; Reynolds & Ogdcn 1979). The velocity gradient of 



Filament 1 is outside the range expected from a simple model of Galactic rotation and therefore may be 
evidence for a dynamical influence on the filament. However, the observed velocity gradient is inconsistent 
with this feature being the edge of an expanding shell, since a shell's projected edges should be at a constant 
radial velocity. Furthermore, for the two cases presented here, we find no obvious correlation between these 
Ha filaments and emission at other wavelengths, including H I 21 cm, Rosat All-Sky Survey X-Ray, and 
IRAS bands, making it difficult to relate them to previously identified H I "worms" and superbubbles. Since 



the Leiden/Dwingcloo H I survey is sensitive to below 5 x 10 cm ( Harmann fc Burton 1997 ), the ionized 
column density of the filament, 2 x 10 19 cm~ 2 , suggests that the gas is fully ionized. Other possibilities 
include a suggestion by Dupree & Raymond (1983) that faint, ionized trails of ionized hydrogen could 
be produced by high-velocity white dwarfs. However, the recombination time discussed above suggests 
unreasonable velocities for the star. 

In summary, we present the discovery of two long, faint Ha filaments at high latitudes from the 
WHAM Ha survey. Their origin is not yet identified, but the existence of such features may help to explain 
processes responsible for the maintenance of the general WIM layer. As additional portions of the Ha 
survey are reduced, new clues about the nature of these structures may be revealed. We plan to follow 
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these observations with [S II] and [N II] observations of this region. Additional observations will provide 
information about the gas's temperature and ionization state, which are useful in narrowing the kinds of 
processes that can be producing these filaments. 

We thank Kurt Jaehnig and Jeff Percival of the University of Wisconsin's Space Astronomy Lab 
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Fig. 1. — A portion of the WHAM Ha Sky Survey. These pseudo-color images show the integrated Ha 
emission in four selected velocity bands. The axes are Galactic longitude and latitude, (o) t>LSR = ~60 to 
+60 km s _1 . (6) v hSR = -26 to -14 km s _1 . (c) v LSR = -6 to +6 km s _1 . (d) v LSR = +14 to +26 km s _1 . 
"Holes" in the images are pointings uncorrected for a bright stellar absorption line. (Preprint Note: A color 
PostScript version of this figure can be found at http://www.astro.wisc.edu/wham/papers.html) 

Fig. 2. — Longitude- Velocity intensity map at b = +25°. Contour levels are in units of 10 -3 R (km s -1 ) -1 . 
Filament 1 crosses b = +25° at I = 225°. 

Fig. 3. — Filament and background spectra at b = +25?5. Intensity is plotted against ulsr- Atmospheric 
lines have been removed. The dashed line is the original emission toward the filament, the dotted lines are 
two background directions, and the solid line is the background-subtracted filament emission (see text). The 
increased noise near —40 km s _1 is a residual of the geocoronal Ha line, whose peak emission is centered 
near that velocity. 
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Tabic 1. Fitted Filament Parameters 
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a This direction is toward the H II region surrounding CMa 
R1/OB1. 
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